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INTRODUCTION 


This  final  report  covers  the  period  August  1,  1977  to 
September  30,  1980.  It  really  isn't  final,  since  work  will  continue 
under  a  new  contract  (DAAG29-80-K-0019)  which  carries  the  title 
"Fundamental  Electronic  Properties  of  Donor-Type  Graphite  Intercalated 
Compounds" 

Since  all  the  significant  results  have  been  published  in  the  open 
literature,  we  present  here  the  highlights  of  those  results.  Attached 
as  appendices  are  preprints  of  papers  which  are  not  in  print  as  of 
this  writing. 


STATEMENT  OF  PROBLEM 

-^Graphite  intercalation  compounds  are  a  large  family  of  quasi- 
two  dimensional  synthetic  metals.  Fundamental  interest  centers  on 
their  unusual  crystal  structures  and  their  effects  on  electronic 
properties.  Practical  applications  are  envisaged  as  lightweight 
conductors  of  electricity. 

'"'^The  goals  of  this  program  are  to  understand  the  electronic 

structure  of  donor- type f compounds ,  with  a  view  toward  obtaining 

detailed  correlations  between  chemical  and  physical  variablles. 

Emphasis  is  placed  on  the  one-electron  spectrum,  transport  and 

optical  properties  and  charge  distributions.  Synthetic  methods  are 

1  • 

developed  to  provide  high-quality  specimens  for  a  variety  of  physical 
measurements.  Band  theoretical  techniques  are  applied  to  develop 
models  for  the  electronic  structure.  Iteration  between  theory  and 
and  experiment  leads  to  more  sophisticated  models./ 


HIGHLIGHTS 


1.  It  had  been  known  for  some  time  that  the  Hall  coefficient 

of  KCg  and  related  materials  was  anomalously  small  at  room 

temperature,  implying  a  free  electron  density  in  excess  of  the  metal 

atom  density  if  one  chose  to  interpret  Ru  with  a  one-carrier  model. 

n 

We  found  from  low  temperature  measurements  that  the  hall  coefficient 
changes  sign,  becoming  p-type  in  the  helium  range.  This  suggested  a 
two-carrier  model  for  the  transport  process,  in  which  electrons  and 
holes  are  nearly  compensated  at  room  temperature.  For  the  reflectance 
spectra,  a  two-carrier  model  was  also  required  to  obtain  a  good  fit  to 
the  data,  as  well  as  to  obtain  agreement  between  optical  parameters 
and  the  dc  conductivity.  This  interpretation  remains  somewhat 
controversial,  in  the  sense  that  the  Japanese  group  under  Prof.  Tanuma 
do  not  find  a  positive  Hall  coefficient  at  low  temperature.  We  simply 
point  out  that  the  resistance  ratios  of  our  samples  are  of  order  200- 
300,  while  Tanuma' s  samples  are  of  order  10. 

2.  With  the  publication  of  band  structures  for  KCg  and  LiCg,  by  the 
Japanese  and  ourselves  respectively,  it  became  apparent  that  the 
different  in-plane  symmetries  of  intercalant  layers  in  these  two 
materials  might  be  playing  an  important  role  in  determining  the  overall 
electronic  properties.  In  KCg,  the  metal  s  electrons  were  found  to 
hybridize  strongly  with  the  carbon  pz  electrons  near  the  Fermi  energy, 
giving  rise  to  two  types  of  carriers  (in  agreement  with  the  data 
mentioned  above).  In  LiCg,  on  the  other  hand,  no  evidence  was  found 
for  s-pz  mixing  near  Ep.  At  the  time,  LiCg  was  a  difficult  material 


experimentally,  due  to  limitations  in  the  available  synthetic 
methods.  We  proceeded  to  develop  simplified  techniques  for  making 
LiCg  and  the  stage  two  compound  "LiC^"/  and  to  supply  samples  of  these 
materials  to  many  groups  for  a  wide  variety  of  experiments. 

3.  We  found  that  the  temperature  dependence  of  resistivity  in  the 

heavy  alkali  compounds  could  be  fit  to  the  sum  of  a  linear  and  a 

quadratic  term,  the  linear  term  representing  usual  electron-phon 

scattering  while  the  quadratic  term  was  interpreted  as  either  electron- 

hole  or  electron-electron  scattering  in  a  complex  Fermi  surface.  This 

interpretation  agreed  with  the  band  structure  work  noted  above,  in 
2 

that  the  T  term  was  much  larger  in  stage  one  materials  than  in  stage 
two. 

4.  In  collaboration  with  DiSalvo  and  Safran  at  Bell  Laboratories,  we 
performed  total  susceptibility  measurements  on  stages  one,  two,  three 
and  four  potassium  samples.  Surprisingly,  we  found  that  the  suscepti¬ 
bility  for  magnetic  fields  along  the  c-axis  increased  with  stage, 
contrary  to  what  would  be  expected  if  the  origin  of  the  susceptibility 
were  simply  Pauli  paramagnetism.  This  effect  was  ultimately  traced  to 
an  unusual  Landau  paramagnetism  which  turns  out  to  be  characteristic 
of  "doped"  two-dimensional  graphite.  In  other  words,  if  in  a  high 
stage  compound  the  transferred  electrons  are  mainly  localized  on 
carbon  layers  neighboring  the  metal  layers,  then  these  carbon  layers 
will  exhibit  an  anomalously  large  orbital  paramagnetism.  The  interior 
layers,  those  far  from  a  metal  layer,  are  more  nearly  like  pure  graphite, 
and  should  ultimately  (at  some  very  high  stage)  exhibit  the  large  orbital 


diamagnetism  associated  with  the  degenerate  semimetallic  nature 
of  pure  graphite.  It  was  realized  that  by  extending  these  measure¬ 
ments  to  more  dilute  materials,  one  might  be  able  to  map  out  the 
c-axis  charge  distribution  in  detail. 

5.  We  performed  a  detailed  band  structure  calculation  for  LiC^ , 
choosing  this  compound  for  its  simple  crystal  structure.  A  Green's 
function  technique  with  non-muf f in-tin  corrections  was  used,  and 
approximate  self-consistency  was  checked  by  comparing  results  from 
atomic  and  fully  ionic  potentials.  This  calculation  is  now  widely 
accepted  by  workers  in  the  field,  because  of  the  many  aspects  in  which 
it  agrees  with  experiment: 

-  The  calculated  density  of  states  at  the  Fermi  energy  N(Ep) 
agrees  with  the  low  temperature  linear  term  in  specific  heat  to  better 
than  one  percent. 

-  The  metallic  reflectivity  spectrum  is  fit  very  well  by  the 
theoretical  plasma  frequency  and  a  reasonable  core  dielectric  constant. 

-  The  anisotropy  in  the  conductivity  at  room  temperature  agrees 
well  with  the  calculated  anisotropy  of  the  plasma  frequency,  implying 
that  the  scattering  time  is  essentially  isotropic. 

-  The  experimental  values  for  Knight  shift  and  Pauli  susceptibility 
are  consistent  with  the  theoretical  value  for  the  contact  density  of 

s  electrons. 

-  The  angle  resolved  photoemission  spectrum  of  LiC^  is  in  excellent 
agreement  with  the  calculated  energy  band  structure. 


6. 


Li 


6.  The  theoretical  techniques  developed  for  LiCg  are  being  applied 
to  the  first-stage  potassium  compound  KCg.  The  energy  band  picture 
is  substantialy  completed  and  the  results  indicate  that  the  bottom 
of  K  4s  band  lies  ^1.5  eV  above  the  Fermi  level.  This  indicates  that 


the  Fermi  surface  will  lack  the  central  isotropic  pocket  of  carriers 
predicted  by  other  studies.  A  series  of  saddle  point  transitions  centered 
at  4.1  eV,  appear  consistent  with  the  low  energy  electron  loss  spectra. 

7.  A  similar  theoretical  study  of  pur  graphite  was  initiated  to  test 
the  technique  used  for  LiCg  and  KCg  on  a  better  known  material.  The 
initial  picture  of  the  energy  bands  is  in  good  agreement  with  those 
obtained  by  other  techniques. 
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ABSTRACT 


Elastic  and  inelastic  neutron  scattering  are  reported  for 
stage  1  and  2  Li-graphite.  The  principal  results  are  as  follows 

1.  No  evidence  was  found  for  the  previously  reported 
stacking  transition  in  LiC^;  the  normal  AaAu  ... 
stacking  persists  down  to  10K. 

2.  The  C-parameter  of  LiC^  increases  rapidly  with 
increasing  T  (1.5%  between  300  and  490K) .  This 
may  be  a  precursor  to  the  Li  order-disorder 
transition  which  occurs  above  500K. 

3.  Upon  heating  above  520K,  stage  2  transforms  into 
a  dilute  stage  1  phase,  as  predicted  by  Safran. 

4.  The  dispersion  of  the  LA  phonon  propagating  along 
c  shows  a  40%  increase  in  sound  velocity  relative 
to  pure  graphite,  as  compared  to  10%  for  RbCg. 

This  confirms  the  existence  of  strong  interlayer 
interactions  in  LiCg . 
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INTRODUCTION 


Several  papers  at  this  conference  deal  with  phase  transitions 

in  various  compounds.  Such  studies  are  important  both  in  terms 

of  expanding  our  knowledge  of  intercalated  graphite'*'  and  as  they 

2 

apply  to  general  issues  in  critical  phenomena.  Bak  has  demonstrated 
the  isomorphisms  between  prototypical  GIC's  and  idealized  model 
lattices  whose  critical  behavior  have  been  predicted  theoretically. 
Experiments  on  GIC's  therefore,  can  serve  as  a  basis  for  testing 
these  theories. 

We  report  here  the  preliminary  results  of  a  neutron  scattering 
study  of  phase  transitions  and  phonon  spectra  in  stages  1  and  2 
Li-graphite.  Large  (''■I  gram)  samples  were  prepared  from  HOPG  by 
liquid  phase  synthesis3  and  sealed  in  quartz  or  stainless  steel 
tubes  along  with  He  exchange  gas.  Since  the  isotope  6Li  (4%  abundant) 
is  a  strong  absorber  for  thermal  neutrons,  we  used  99+%  enriched 
^Li.  Neutron  scattering  experiments  were  performed  on  a  triple 
axis  spectrometer  at  the  Siloe  reactor  in  Grenoble.  Pyrolitic 
graphite  was  used  as  monochromator,  analyzer  and  filter.  A  wave¬ 
length  of  2.4  A  was  used  for  the  elastic  diffraction  studies. 

The  sample  was  placed  in  either  a  He  cryrostat  or  an  oven,  with 
c-axis  in  the  scattering  plane. 
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SEARCH  FOR  STACKING  TRANSITION 

At  300K  the  Li  atoms  in  LiC^  form  a  triangular  lattice 

commensurate  with  the  hexagonal  array  of  carbon  atoms.  The  unit 

cell  is  hexagonal  with  a=aolfT,  where  aQ  is  the  graphite  in-plane 

parameter.  All  Li  layers  are  equivalent,  as  are  all  C  layers, 

so  the  alternate  stacking  in  stage  1  can  be  written  AaAn . . . , 

4  5 

A  denoting  carbon  and  a  denoting  Li  layers.  Kambe  et .  al. 
recently  reported  a  transition  to  a  more  complex  stacking  upon 
cooling  below  200k.  They  interpreted  their  electron  diffracto- 
grams  for  T<200K  in  terms  of  AaABAy .  .  .  .  stacking,  based  on  the 
disappearance  of  (h,Q,0)  diffraction  spots,  the  only  remaining 
spots  being  (h  0,0),  i.e.  those  due  to  the  graphite  in-plane 
structure . 

With  a  view  toward  studying  this  transition  in  more  detail, 
we  recorded  (0,6, si),  (h,0,0)  and  (1,0,  l)  diffraction  scans  at 
10,  200  and  300k.  To  our  surprise,  no  transition  was  observed. 

A  sampling  of  the  data  is  shown  in  Fig.  1.  The  top  panel  shows 
a  (1,0, f)  scan  at  10K.  The  Bragg  peaks  at  d=Q,  1  are  characteristic 
of  the  hexagonal  cell,  i.e.  AaAot...  stacking.  If  a  transition  to 
AaABAy  had  occurred,  the  unit  cell  would  now  be  rhombohedral , 
leading  to  extintion  of  the  1=0,  1  peaks  and  the  appearance  of  new 
peaks  at  £=1/3,  2/3  (both  occurring  due  to  the  powder  averaging 
of  in-plane  scattering  in  HOPG) .  The  middle  and  lower  panels  show 
(h,0,0)  scans  at  200K  and  10K.  Here,  the  fact  that  the  (1,0,0) 
and  (2,0,0)  intensities  do  not  change  upon  cooling  also  signifies 
the  absence  of  a  stacking  transition. 
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For  butt  scans  the  300K  data  are  identical  to  the  10K  results,  and 

two  samples  behaved  identically.  We  conclude  that  the  ground  state 

structure  is  the  same  as  the  300K  one  first  reported  by  Guerard  and 
4 

Herold,  and  that  some  other  phenomenon  must  be  invoked  to  explain 
the  electron  diffraction  results^ 

It  must  be  pointed  out,  however,  that  despite  the  discrepancy 
between  electron  and  neutron  diffraction  results,  weak  anomalies 
have  been  observed  in  specific  heal?  and  basal  plane  conductivity 3at 
temperatures  corresponding  to  the  change  in  electron  diffraction  pattern. 
On  the  other  hand,  no  effect  was  found  in  differential  scanning 
calorimetry^ 

HIGH-TEMPERATURE  RESULTS 

Diffraction  scans  were  recorded  for  a  LiC^  sample  sealed  in 
quartz,  as  a  function  of  increasing  temperature.  Unfortunately  the 
container  failed  before  reaching  the  order-disorder  transition, 
evidently  due  to  a  trace  of  free  Li  which  attacked  the  quartz. 

However,  the  elastic  scans  revealed  a  very  large  thermal  expansion 
in  the  c-direction,  shown  in  Fig.  2.  The  expansion  is  nonlinear, 
and  the  rate  increases  with  increasing  T.  Over  the  same  interval 

-6 

t  tv ■  a-parameter  is  constant  within  experimental  error  (Aa/a  5X10/K  )  . 

.'h ;  s  suggests  anomalously  large  out-of-plane  thermal  motion  of  Li 
as  the  order-disorder  transition  is  approached. 

A  second  sample  was  packaged  in  a  stainless  steel  tube  thinned 
to  0.007  inches  to  minimize  the  background  scattering.  This  sample 
also  suffered  from  an  unplanned  event  during  heating,  namely  the 
failure  of  a  thermocouple  which  led  to  a  temperature  excursion  above 
650K. 
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Upon  cooling  to  room  temperature,  it  was  found  that  the  sample  had 

irreversibly  transformed  to  a  mixture  of  stage  2  (presumably  LiC^) 

and  lithium  carbide  Li2C2*  After  repairs,  the  stage  2  phase  was 

studied  as  a  function  of  increasing  T.  Figure  3  shows  the  temperature 

variation  of  the  integrated  (0,0,1)  peak  of  "LiC^"  (this  corresponds 

to  2=0.52  in  units  of  the  LiCg  cell).  The  transition  is  clearly  a 

higher  order  phenomenon  with  large  hysteresis.  The  remaining 

intensity  at  150-200C  is  due  to  a  new  expanded  stage  2  phase 

(LiClg?)  .  By  250C  the  only  remaining  (0,0,2.)  intensity  is  indexable 

as  stage  1,  but  with  c  parameter  ^1%  smaller  than  that  of  Lie, 

o 

at  that  temperature .  Scans  along  (h,0,0)  at  250c  show  that  the 
Li  layer  is  now  disordered.  We  conclude  that  the  net  effect  over  the 
range  100C  -  250C  is  the  transformation  of  ordered,  close-packed 
stage  2  to  dilute  stage  1  with  liquid-like  Li  layers,  as  predicted 
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by  Safran.  A  dilute  stage  2  intermediate  phase  seems  to  occur  as 
the  stage  1  grows  at  the  expense  of  "dense"  stage  2.  This  effect 
is  completely  reversible.  We  plan  to  repeat  this  experiment  on 
a  well-characterized  LiC^  sample. 

PHONONS 

Inelastic  neutron  scattering  of  LiC^  was  performed  at  200K 
in  a  geometry  which  reveals  only  longitudinal  modes  propagating 
in  the  c-direction.  Figure  4  shows  the  loss  spectra  for  two  q-values. 
Two  loss  peaks  are  observed,  both  of  which  increase  in  energy  with 
increasing  q.  The  lower  energy  peak  decreases  in  intensity  with 
increasing  q,  and  is  no  longer  observed  beyond  q=0 . 3  ( 2ir/c)  . 


f 

Figure  5  shows  the  dispersion  of  these  two  modes.  The  upper  mode 
extrapolates  to  zero  energy  at  q=0,  so  we  assign  it  to  the  LA 
phonon.  On  the  other  hand,  the  lower  mode  tends  toward  a  value 
1-2  meV  (8-16K)  at  q=0,  so  we  tentatively  assign  this  branch  to  an 
LO  mode . 

The  LA  mode  can  be  compared  directly  to  the  corresponding 
8  9 

branches  in  graphite  and  RbCg.  From  the  initial  slopes  we  find 
that  the  c-axis  sound  velocity  of  LiCg  is  enhanced  40%  relative 
to  graphite,  whereas  RbCg  is  enhanced  only  10%.  This  in  turn  implies 
a  relatively  large  interlayer  interaction  in  LiCg ,  which  is  consistent 
with  many  independent  observations: 

1)  The  observed  layer  spacing  is  only  83%  of  the  hard-sphere 
estimate,  suggesting  strong  covalency? 

2)  Theoretical  charge  density  indeed  reveals  a  pileup  between 
11 

C  and  Li  layers. 

3)  The  conductivity  anisotropy  oa/oc  is  only  10  at  300K? 

6  3 

4)  The  lattice  specific  heat  goes  as  T  (i.e.  three-dimensional) 
up  to  50K,  whereas  in  graphite  the  T'*  region  extends  only 
to  MK. 

The  putative  LO  mode  is  harder  to  explain.  No  similar  mode 
was  found  in  RbCg,  so  we  are  led  to  suspect  that  the  small  mass  of 
Li  is  implicated.  One  possibility  would  be  a  ferroelectric  instability 
whereby  at  low  T  the  Li  atoms  are  located  in  two  local  potential 
minima  straddling  the  mid-plane. 


The  barrier  between  these  two  wells  would  correspond  to  the 
q=0  "gap"  of  the  LO  branch,  i.e.  ^10K.  Above  10K,  thermal  motion 
would  produce  a  structure  in  which,  on  average,  the  Li  layer  is 
equidistant  between  flanking  carbons,  as  observed.  This  hypothesis 
is  subject  to  verification  by  looking  for  a  displacive  transition  of 
Li  layers  below  10K,  and  by  determining  the  temperature  dependence 
of  the  LO  mode. 

A  second  possibility  is  that  the  mode  results  simply  from  '.he 

Li-carbon  interaction,  which,  among  other  things,  could  mix  some  of 

8 

the  low-lying  modes  which  occur  in  graphite.  Calculations  show  that 
the  TO  and  LA  modes  in  pure  graphite  cross  each  other  halfway  from 
q=0  to  the  zone  boundary.  It  is  conceivable  that  in  LiCg  these 
graphitic  branches  produce  a  low-lying  mode  with  a  small  gap  at 
q=0.  This  hypothesis  can  be  tested  by  a  lattice  dynamics  calculation 
for  LiCg 

We  acknowledge  fruitful  discussions  with  D.  Guerard,  J.Y.  Henry 
and  S.  A.  Solin. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 
Figure  3 


Figure  4 


Figure  5 


Elastic  scans  for  LiC^  at  various  temperatures  and  in 
various  directions,  as  indicated.  Indexing  based 
on  the  Guerard  cell  (Ref.  4),  with  a=4.30A,  c=3.7lA. 

No  stacking  transition  is  observed  down  to  10K 
(see  text) . 

C-axis  thermal  expansion  of  LiC^ ,  as  measured  by  the 
temperature  variation  of  the  (0,0,1)  Bragg  peak. 

Temperature  dependence  of  the  integrated  (0,0,1)  peak  of 
stage  2,  LiC^*  The  sample  eventually  trangforms 
entirely  into  a  dilute  stage  1  phage  at  250°C.  The 
remaining  stage  2  intensity  at  200  C  is  due  co  a 
dilute,  slightly  expanded  stage  2  intermediate  phase. 

Inelastic  neutron  energy  loss  spectra  of  LiC,  at 
200K,  recorded  for  q=0.1,  0.2  (the  zone  boundary  is 
at  q=0.5).  These  are  longitudinal  modes  propagating 
along  c. 

Dispersion  of  the  two  modes  in  Fig.  4.  The  upper  branch 
is  the  LA  mode,  while  the  lower  branch  is  tentatively 
assigned  to  an  LO  mode  (see  text) . 


Neutron  counts  x  10 


#/.«/  />*}  I 


Scans  along  [  hOO]  Scan  along  [1011 


Neutron  ^counts  (*0.01) 


1 


Energy  (meV) 


(Cj'j  £f- 


Elsevier  Sequoia  S.A  .  P.O.  Box  851.  1001  Lausanne  1.  Switaerland 

SYNTHETIC  METALS 


XfioiTOH 


j  CORRECTION 


mastercopy[ 


'AUTHOR 

QUERIES 


Only  typographical  correction  will  bs  accepted  at  thii  itagei 


EZ] 


Vol. 


IN  ALL  COHHESPONDENCE 
CONCERNING  THIS 
PAPEH  RtfER  TO 


SYM 


Synthetic  Metals, 

©  Elsevier  Sequoia  S.A.,  Lausanne  —  Printed  in  the  Netherlands 


A  CALCULATION  OF  THE  ELECTRONIC  BAND  STRUCTURE  OF  KC8 

i  D.  P.  DIVICENZO*,  N.  A.  W.  HOLZWARTH  and  S.  RABII 

Moore  School  of  Electrical  Engineering,  University  of  Pennsylvania,  Philadelphia,  PA 
19104  (U.S.A.) 

(Received  May  19, 1980) 


Summary 

Results  of  an  ab  initio  band  structure  calculation  for  first  stage  potas¬ 
sium-graphite,  KCg,  are  presented  along  several  high  symmetry  directions  in 
the  Brillouin  zone.  In  the  first  step,  we  find  the  KKR  band  structure  using 
the  muffin  tin  approximation.  In  the  second  step  we  introduce  the  correc¬ 
tions  to  the  muffin  tin  potential,  leading  to  the  final  bands. 

Our  band  structure,  excluding  the  bands  originating  from  the  3p  and  4s 
levels  of  potassium,  is  in  good  agreement  with  the  predictions  of  the  folded 
band  model  of  two  dimensional  graphite.  The  bottom  of  the  K  4s  band  lies 
about  1.5  eV  above  the  Fermi  level.  Preliminary  indications  are  that  our 
Fermi  surface  will  lack  the  central,  isotropic  pocket  of  carriers  predicted 
earlier,  and  will  be  more  in  agreement  with  the  folded  bands.  We  have  also 
studied  the  lowest  energy  critical  point  transitions  and  find  a  series  of  saddle 
point  transitions  centered  at  4.1  eV,  roughly  consistent  with  the  experimen¬ 
tal  low  energy  electron  loss  spectra. 


1.  Lotroduction 

-  First  stage  potassium-graphite,  KC8,  continues  to  be  a  system  under  ac¬ 
tive  experimental  and  theoretical  investigation.  In  the  last  year  alone,  mea¬ 
surements  of  electron  energy  loss  [1],  angle  integrated  [2]  and  angle  re¬ 
solved  [3]  photoemission,  de  Haas-van  Alphen  periods  [4],  and  electrical 
conductivity  and  Hall  coefficient  [5]  have  been  performed.  On  the  theory 
side,  calculations  of  the  electronic  band  structure  and  Fermi  surface  proper¬ 
ties  [6]  and  a  KC8  phonon  spectrum  [7]  have  been  reported.  The  proximity 
of  |>otassium  states  to  the  Fermi  level  of  KC8  makes  an  accurate  band  cal- 
culation  difficult;  we  hope  that  our  electronic  structure  work,  based  on  a 
mathematically  accurate  ab  initio  technique,  can  provide  a  complete  and  ac¬ 
curate  theoretical  picture. 
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This  paper  serves  as  a  brief  progress  report  on  our  KCg  band  calcula¬ 
tion.  Section  2  presents  KKR  muffin  tin  bands  throughout  the  Brillouin 
zone.  Section  3  develops  the  method  which  we  use  to  correct  the  muffin  tin 
approximation,  discusses  some  of  the  difficulties  which  we  encountered,  and 
gives  the  final  bands  obtained  to  date.  Section  4  relates  our  results  to  some 
of  the  experiments  mentioned  above,  considers  this  calculation  in  relation  to 
previous  ones,  and  speculates  on  work  to  be  done  in  the  near  future. 


2.  Muffin  tin  energy  bands 

Our  method  of  band  structure  calculation  has  been  described  previously 
[8,9],  The  plan  of  our  work  is  as  follows:  A  crystal  potential  is  formed  by 
the  linear  superposition  of  atomic  charge  densities  for  K+1  and  CT1/8.  In  the 
first  step  we  use  a  muffin  tin  approximation  to  this  potential;  a  constant  out¬ 
side  the  atomic  spheres,  the  actual  value  inside.  YVe  calculate  the  scattering 
matrix  kr lml-m  -and  the  KKR  structure  matrix  MTr  then  solve  the  secular 
equation: 

'  det  jfe-1  wm  •(£)& rr-  +  >(k  ,£)J  =  0  (1) 

P  [vor  wavevectors  throughout  the  Brillouin  zone.  The  result  of  this  calculation 
is  shown  in  Fig.  1. 

Although  the  results  of  Fig.  1  are  only  approximate,  they  contain  some 
physically  useful  information.  First,  we  have  discovered  that  for  low  energies 
(below  about  —0.3  Ry)  the  wave  functions  have  relatively  unmixed  potas¬ 
sium  or  carbon  character,  so  that  the  set  of  dispersionless  bands  at  —1.39  Ry 
is  clearly  identifies  as  K  3p,  the  pair  of  bands  with  minimum  at  —0.48  Ry. 

(T  j*)  and  —0.37  Ry  (r1~)  are  derived  from  K  4s  states,  while  the  rest  of  the 
bands  come  from  the  states  of  pure  graphite.  Moreover,  these  graphite  levels 
-  agree  in  marui  respects  with  those  predicted  by  the  folded  band  model.  Fig¬ 
ure  2  shows  the  energy  bands  of  two  dimensional  graphite  as  calculated  by 
Nagayoshi  et  al.  [10]  folded  into  the  Brillouin  zone  of  KC8.  A  comparison 
of  Figs.  1  and  2  shows  that  the  graphite  a  bands  (shown  solid)  and  most  of 
the  graphite  rr  bands  (dashed)  survive  with  very  little  distortion  in  potas¬ 
sium-graphite.  Also,  all  of  the  predicted  band  degeneracies  shown  in  paren- 
^  theses  in  Fig.  2  are  correct  for  the  KC8  bands. 

Still,  there  are  certain  aspects  of  the  KKR.  muffin  tin  calculation  which 
change  significantly  as  a  result  of  the  muffin  tin  corrections.  When  compared 
with  the  actual  crystal  potential,  the  muffin  tin  potential  is  too  attractive  in 
the  intercalant  plane  and  too  repulsive  in  the  carbon  plane.  This  means  that 
the  relative  positions  of  the  graphite  n  and  a  bands,  as  well  as  that  of  the 
potassium  and  the  graphite  bands,  are  not  given  properly  by  this  calculation. 
The  potassium  4s  level,  for  example,  comes  roughly  two  eV  below  the  graph¬ 
ite  Fermi  level  at  —4.5  eV,  in  violent  conflict  with  experiment.  Thus  it  is 
pointless  to  attempt  any  further  analysis  —  Fermi  surface  effects,  optical 
properties,  etc.  —  of  our  muffin  tin  bands. 
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Fit'.  1-  The  KKR  muffin  tin  bands  of  potassium  graphite.  Bands  derived  from  graphite  7 T 
states  are  dashed ;  solid  bands  are  graphite  <7  and  potassium  states. 


3.  Final  KC8  bands 

The  deviation  of  the  full  crystal  potential  from  the  muffin  tin  value  in 
the  interstitial  region,  A  V(r ),  may  be  taken  into  account  by  calculating  its 
matrix  element  between  suitably  normalized  muffin  tin  wave  functions!^. 

‘  <*/):  _  ° 
A 


We  now  find  the  final  energy  eigenvalues  of  KC8  by  solving  a  second 
secular  matrix: 

de;  [Ann-(^)  +  (£nKKR(fc')-£)5rol'l  =  0  (3) 

...  yrheroJLK,KR  are  all  the  energy  eigenvalues  calculated  in  the  last  Section 
within  a  range  of  3Ry. 

The  evaluation  of  the  matrix  elements  A nn-(ft)  requires  some  discus¬ 
sion.  For  most  high  symmetry  wavevectors  the  integral  (2)  need  only  be 
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Fig.  2.  Bands  of  two  dimensional  graphite  folded  into  the  Brillouin  zone  of  KC8.  n  Bands 
are  dashed.  O  bands  solid.  Numbers  in  parentheses  indicate  the  band  degeneracies.  The 
insert  shows  the  relationship  between  the  graphite  and  KCg  Brillouin  zones. 


evaluated  in  |  of  the  unit  cell,  the  irreducible  sector.  Also,  the  integral  is 
only  non-zero  in  the  interstitial  region,  which  makes  up  only  43%  of  the 
crystal  volume.  Howe/er,  the  irregularity  of  this  integration  region  (see  Fig. 
3)  precludes  the  use  of  any  of  the  standard  numerical  integration  methods. 
We  have  developed  a  grid  which  approx  mates  the  integral  (2)  using  Gaussian 
quadrature  [11]  in  the  three  coordinate  directions  and  which  avoids  the  dis- 
-continuities  caused  by  the  irregular  geometry.  After  considerable  testing  an 
integration  grid  containing  1675  points  in  the  interstitial  portion  of  the  ir¬ 
reducible  sector  was  used.  Our  previously  published  work  [9]  used  a  483 
point  grid,  which  was  not  completely  converged.  Thus,  the  results  reported 
were  not  absolutely  correct.  In  our  most  recent  work,  numerical  errors  in 
A „„•(£* )  are  no  greater  than  0.03  Ry  [12] . 

W’e  have  finished  the  solution  of  eqn.  (3)  for  several  wavevectors.  Our 
results  are  shown  in  Fig.  4.  Again,  we  show  the  graphite  a  and  potassium 
bands  solid,  the  -  bands  dashed.  The  KC8  band  shapes  still  bear  a  close 
resemblance  to  the  folded  graphite  bands  of  Nagayoshi  et  al.  in  Fig.  2.  The 
relative  positions  of  the  graphite  a  and  n  bands  now  agree  very  closely  with 
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Fit'  3.  KCg  irreducible  sector  showing  inscribed  muffin  tins  of  potassium  and  carbon. 


the  bands  of  two  dimensional  graphite  folded  into  the  KC8  Brillouin  zone. 
The  potassium  3p  levels  appear  as  the  flat  set  of  bands  at  —1.39  Ry.  The  pair 
of  bands  with  minima  at  —0.26  Ry  (r^)  and  —0.15  Ry  (r1~)  are  derived 
from  the  K  4s  levels.  As  expected,  the  potassium  bands  have  shifted  dramat¬ 
ically  from  their  positions  in  the  muffin  tin  calculation. 


4.  Discussion 

Although  we  have  completed  our  calculation  for  only  a  few  U  points 
in  the  Brillouin  zone,  we  may  draw  several  reasonable  conclusions  from 
them.  The  KC8  Fermi  level  will  fall  at  about  —0.4  Ry  on  Fig.  4.  The  energy 
bands  at  this  energy  and  below  have  very  similar  shapes  to  the  folded  graph¬ 
ite  bands,  although  most  of  them  are  split  by  interactions  with  potassium 
states.  We  conclude  that  the  density  of  states  will  be  very  similar  to  graph¬ 
ite’s,  and  the  Fermi  surface  will  agree  with  that  derived  from  the  rigid  band 
analysis  performed  by  Holzwarth  et  cl.  [3] .  This  surface  consists  of  rounded 
trigonal  prisms  centered  at  K  and  equivalent  wavevectors  in  the  Brillouin 
zone.  This  differs  from  the  Fermi  surface  calculated  by  Inoshita  et  al.  [6]  in 
"  ""  that  it  lacks  the  roughly  isotropic  pocietol  carriers  centered  at  the  T  point; 
detuned  calculations  will  be  required^to  see  if  existing  experimental  results 
car.  oe  explained  by  our  Fermi  surface. 

Table  1  gives  the  energies  of  the  lowest  energy  critical  point  interband 
transitions  at  f,  which  are  not  critically  sensitive  to  the  exact  position  of  the 
Fermi  level.  All  the  transitions  shown  derived  from  the  M-  point  transition  in 
the  -  bands  of  graphite  reported  at  4.5  eV  [1],  which  is  broadened  and 
lowered  in  energy  by  the  presence  of  the  intercalant.  Hwang  et  al.  [1]  have 
studied  KC„  using  electron  energy  loss.  They  report  a  broad  peak  in  e2(£) 
centered  at  3.8  eV,  in  rough  agreement  with  our  results. 

Future  calculation  should  allow  us  to  perform  a  more  detailed  analysis 
of  Hwang’s  data,  and  to  study  other  experimental  results.  As  an  example, 
Oelhufen  ct  al.  (2  j  have  interpreted  their  photoemission  data  in  terms  of  a 
partially  filled  potassium  4s  band,  in  disagreement  with  our  results.  We  hope 
to  proride  a  new  interpretation  of  their  data. 
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Fig.  4.  Preliminary  results  of  final  corrected  band  structure  of  potassium  graphite.  Jr 
Bands  are  dashed,  o  bands  and  potassium  levels  are  solid. 
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TABLE  1 

Optical  selection  rules  and  principle  low  energy  critical 
transitions  for  the  T  point  in  KCg  ^ 


II 

II 

1 

1+  *♦  3— 

3+  **  2— 

2+ 

l- 

II 

1 

n 

1+  **  4  — 

3+  ++  4— 

44 

4-k 

2- 

A 

II 

u 

1+  ♦*  2— 

2+-  ♦*  4— 

4  + 

«-> 

1— 

II 

u 

A 

3+  -»  1- 

2+  **  3— 

4  + 

4-> 

3- 

Transition  at  T 

A E  (eV) 

4 - *  2+ 

3.4 

4 —  !♦ 

3.5 

3 - *  2+ 

3.5 

4 +  -*•  1— 

3.8 

3  +  1  +- 

3.9 

3  r  -*•  1  — 

4.r 

4+  -*  1  — 

4.2 

3+  -►  2— 

4.3 

4*  -*  2— 

4.5 

3 - »  2+ 

4.7 

4-  -*  2  + 

4.8 
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Summary 

We  present  a  progress  report  of  our  ab  initio  energy  band  structure  cal¬ 
culation  of  3-dimensional  graphite.  The  motivation  for  this  work  derives 
from  2  principle  objectives:  (1)  to  provide  a  basis  of  comparison  for  our 
work  on  graphite  intercalation  compounds  and,  (2)  further  to  test  some 
computational  aspects  of  our  formalism.  These  considerations  are  discussed 
further  in  Section  1.  In  Section  2  we  describe  the  graphite  structure  and  our 
crystal  potential.  The  formalism  which  has  been  presented  in  detail  else¬ 
where  [1]  is  briefly  reviewed  in  Section  3.  In  Section  4  we  present  the  re¬ 
sults  of  our  KKR  calculation  for  a  muffin  tin  potential.  The  final  bands  we 
have  obtained  so  far  are  presented  in  Section  5  and  we  close  with  a  compari¬ 
son  of  our  results  with  optical  experiments  and  other  theoretical  studies. 


1.  Introduction 

There  have  been  many  band  structure  calculations  of  graphite  to  date. 

The  band  structure  of  2-dimensional  graphite  has  been  presented  in  refs.  2  - 
6.  The  two  dimensional  model  of  graphite  necessarily  ignores  interlayer 
interactions  which  introduce  band  splittings  of  roughly  1  eV  and  introduce 
,  important  structure  near  the  Fermi  level. 

The  previous  three  dimensional  calculations  of  the  graphite  band  struc¬ 
ture  roughly  fall  in  two  categories.  First  there  are  the  band  structures  which 
are  fitted  in  detail  to  experimental  results  over  a  small  energy  range.  These 
provide  accurate  information  about  the  Fermi  surface  and  are  useful  for  cor¬ 
relating  transport  related  measurements  [5,  7]  and  explaining  optical  struc¬ 
ture  at  moderate  energies  [8,  9] .  The  second  category  of  three  dimensional 
calculations  cover  a  larger  energy  range  and  employ  methods  ranging  LCAO<^- 
[10, 11] ,  pseudo-potentials  [12]  and  original  reformulations  of  combined 
OPW  and  tight-binding  techniques  [13].  In  general,  the  accuracy  of  these 
calculations  is  not  sufficient  to  expect  an  adequate  representation  of  the 
Fermi  surface  of  graphite.  However,  in  ref.  13,  a  good  representation  ob- 
tained  by  tuning  adjustable  parameters  of  the  model  potential.  Higher  energy 
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optical  properties  are  also  better  predicted.  Reference  6  provides  a  good 
comparative  review. 

The  work  was  initially  undertaken  to  test  the  applicability  of  our 
procedure  to  graphite  intercalation  compounds.  Because  graphite  represents 
a  structurally  limiting  case,  i.e.,  highest  geometric  anisotropy  of  the  crystal¬ 
line  graphite  compounds,  and  given  the  vast  body  of  literature  on  graphite, 
it  was  felt  that  it  would  provide  the  best  test  of  our  method. 


2.  Structure  and  potential 

Graphite  has  the  Bernal  structure  [14]  (Fig.  1)  with  four  atoms  per 
unit  cell  ,  two  each  of  the  inequivalent  A  and  B  carbons.  The  space  group  of 
graphite  is  P63/mmc  —  the  subscript  clearly  displaying  the  non-sym- 
morphicity.  There  are  twenty  four  operations  in  the  point  group,  which 
means  that  most  of  the  calculations  can  be  carried  out/bbtonsidering  only 
l/24th  of  the  volume  of  the  unit  cell  in  an  irreducible  sector.  The  high 
anisotropy  of  the  crystal  structure  naturally  leads  to  a  high  anisotropy  of 
the  potential.  This  fact  rules  out  any  calculations  based  on  a  spherically  sym¬ 
metric  or  a  purely  muffin  tin  potential. 

Our  potential  is  computed  in  the  following  manner.  First,  a  self-con¬ 
sistent  atomic  charge  density  for  carbon  is  obtained  using  the  Hartree-Fock- 
Slater  formalism  [15] .  This  incorporates  the  Slater  Xa  [16]  approach  to  ac¬ 
count  for  exchange/correlation  effects.  A  linear  superposition  of  atomic 
charge  densities  from  each  site  is  used  to  compute  the  crystal  potential  at 
any  given  point.  The  same  Xa  parameter  [17-^  0.759,  is  used  in  both  the 
atomic  and  crystal  potential  calculations. 

The  resulting  crystal  potential  is  shown  along  three  directions  in  Fig. 

2.  At  a  distance  half  the  nearest  neighbor  spacing,  the  potential  variation  in 
.he  three  directions  shown  is  greater  tnan  1  Rydberg.  Clearly  a  muffin  tin  ap¬ 
proximation  will  be  very  poor,  so  any  use  of  such  an  approximation  must  be 
intermediate. 


c  6  696 A 


Fig  1 .  Graphite  crystal  structure.  •,  A  carbons;  O,  B  carbons. 
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The  muffin  tin  potential  that  we  use  in  the  first  step  of  our  calculation 
is  constructed  from  the  actual  potential  in  the  following  manner.  Within 
touching  spheres  about  each  site  (muffin  tin  radius)  the  actual  potential  is 
used.  In  the  interstitial  region,  the  potential  is  taken  to  be  the  volume  aver¬ 
aged  potential  in  that  region  (muffin  tin  zero).  For  graphite  the  muffin  tin 
zero  is  —1.12  Ryd.  The  muffin  tin  radius  is  1.84  a.u.  and  results  in  a  13% 
fractional  fill  of  the  unit  cell  by  the  muffin  tin  spheres. 


3.  Formalism 


Our  band  structure  calculation  is  a  modified  version  of  Painter’s  Dis¬ 
crete  Variational  KKR  method  [IF].  In  the  conventional  KKR  formalism 
only  the  spherically  symmetric  component  of  the  potential  is  retained.  We 
generalize  this  method  to  include  higher  order  terms  in  the  angular  momen¬ 
tum  representation  of  the  potential.  This  gives  rise  to  a  system  of  angular 
momentum  coupled  differential  equations  within  the  muffin  tin  spheres 
which  we  solve  for  a  reaction  matrix  KTmfm'. 

The  reaction  matrix  is,  in  general,  a  complicated  function  of  energy 
which  we  compute  on  a  coarse  energy  grid.  For  graphite,  we  found  that  the 
reaction  matrix  components  are  very  smooth  functions  of  energy  over  the 
energy  range  considered,  in  contrast  to  previous  calculations  where  singular¬ 
ities  were  present.  (1, 19] .  This  smoothness  makes  possible  rapid  and  very 
accurate  numerical  interpolation  and  extrapolation. 
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Our  KK.il  equation  takes  the  form; 


+  A°°  ■ 

/ltml  m 


\ 


■t.w'r  =0 


Uni  rn 

where  K ~1  is  the  reaction  mateix,  A  is  the  structure  matrix  whose  elements 
are  evaluated  with  Kwald  techniques  and  IV  arc  the  coefficients  in  our  KKR 
wave  function,  a  denotes  the  atomic  site.  The  details  have  been  presented  in 
ref.  1. 

To  solve  the  KKR  equation  we  write  it  as  an  eigenvalue  equation: 
[A'-1,?,,  ■& aa'  *  AZfrn  0  '  =  ME)W'"> 


The  eigenvalues  A(E),  are  functions  of  the  energy  because  the  KKR  matrix 
is  energy  dependent.  The  eigen  energies  which  are  physically  meaningful  cor¬ 
respond  to  the  zeroes  of  A(E). 

The  accuracy  of  the  KKR  wave  functions  is  limited  by  the  spherical 
harmonic  l  truncation.  Experience  shows,  however,  that  in  many  cases  an 
s-p  expansion  is  quite  accurate  and  in  almost  all  situations  an  expansion  to 
d  waves  is  extremely  good.  For  graphite  we  have  included  components  up  to 
1-2. 

In  graphite,  all  irreducible  representations  of  groups  of  k  vectors  have 
explicit  representations  in  terms  of  the  atomic  centered  spherical  harmonics 
to  l  =  2,  making  this  truncation  physically  sufficient  as  well. 

In  a  short  test  of  the  effect  of  the  d-wave  components,  we  recomputed 
the  reaction  matrix  and  KKR  eigen  energies  excluding  the  1=2  terms.  We 
found  that  most  of  the  bands  shifted  slightly  but  the  doubly  degenerate  a 
levels  shifted  more  dramatically,  reducing  the  separation  between  the  oc¬ 
cupied  and  unoccupied  cr  bands  by  more  than  2  eV. 


4.  Muffin  tin  bands 

iiy  solving  the  KKR.  equation  for  many  fe-vectors  along  the  edges  of  an 
irreducible  sector  of  the  1st  Briliouin  zone  (Fig.  3),  we  obtained  bands 
snown  in  Fig.  4.  Already  many  symmetry  dependent  features  of  the  final 
bands  are  visible  such  as  the  touching  r.  bands  at  K  and  the  shape  of  the 
lower  a  and  rr  bands  in  the  kz  =  0  plane  which  are  characteristic  of  graphite. 


Kij!  3.  Hriltouin  zone. 
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Kig.  4.  Muffin  tin  bands  of  graphite  (tt  bands  are  shown  dashed). 


The  splitting  of  bands  that  are  degenerate  in  single  layer  calculations  is 
small  but  quite  apparent.  The  muffin  tin  potential  exaggerates  this  feature 
slightly  and  the  splitting  is  reduced  somewhat  in  the  final  bands.  This  can  be 
explained  as  follows:  Recalling  the  potential  (Fig.  4)  one  sees  that  the  muf- 
fm  tin  zero  lies  between  the  potential  extrema  of  the  directions  parallel  and 
perpendicular  to  the  layers.  Thus,  in  general,  the  c  axis  muffin  tin  potential 
is  more  attractive  than  that  of  the  actual  crystal  and,  conversely,  the  in-plane 
muffin  tin  potential  is  more  repulsive  on  average.  This  means  that  the  inter¬ 
layer  interactions  (principally  between  p,  orbitals)  —  which  are  responsible 
for  the  splittings  —  are  overestimated  by  the  muffin  tin  potential.  So  one  can 
expect  the  overall  shapes  of  the  bands  perpendicular  to  the  layer  planes  to 
be  retained  in  the  final  results,  but  the  dispersions  along  A,  U  and  P  to 
change  more  drastically.  This  appears  to  be  borne  out  by  our  results  so  far. 
The  occupied  a  muffin  tin  band  width  at  T  is  16  eV  compared  with  its  final 
value  of  14.7  cV  —  a  change  of  less  than  10%  —  showing  the  relative  stiffness 
of  the  overall  band  shapes.  At  the  K-pt  however,  the  3.2  eV  splittling  of  the 
muffin  tin  n  bands  is  reduced  to  about  1.6  eV  in  the  final  bands  —  the  dis¬ 
persion  along  P  reduced  by  a  factor  of  two. 
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5.  Final  n  m -muffin  tin  bands 


To  compute  the  final  bands  and  eigenfunctions  we  first  express  the 
wavefunctions  of  a  given  symmetry  as  a  linear  combination  of  the  muffin 
tin  wavefunctions  of  the  same  symmetry.  By  diagonalizing  the  non-muffin 
tin  hamiltonian  —  which  is  just  the  potential  difference  Vcr(r>  —  Vmt(r>,  on 
this  basis  we  arrive  at  a  set  of  eigenvalues  and  eigenfunctions  of  the  one 
electron  problem  for  our  crystal.  This  is  we  solve: 

{Enbnn-  +  Aon']  Cn'  =  eCn 

where  C>v  are  the  coefficients  of  the  new  wavefunction, 

A(1„-  =  S'l'niVr.'  —  Vmt)^„'dr3 

£„„■  is  the  eigen  energy  of  the  corresponding  muffin  tin  wavefunction  and 
E  is  the  eigenvalue  (energy)  of  the  new  wavefunction. 

To  compute  the  matrix  elements  A„„'  in  this  last  step,  the  wavefunc¬ 
tions  must  be  suitably  normalized.  Our  normalization  procedure  is  indirect 
but  very  efficient  and  is  described  in  refs.  20  and  21.  The  final  bands  have 
been  completed  at  F,  M  and  at  points  along  A,  2  and  T.  The  most  obvious 
effect  of  the  non-muffin  tin  corrections  is  a  general  downward  shift  of  ap¬ 
proximately  0.3  Kyd  for  the  entire  energy  bands.  Some  of  the  bands,  how¬ 
ever,  behave  very  differently.  In  particular,  the  sigma  (ri+)  band  moves  up 
from  the  Fermi  level  into  the  higher  conduction  bands  (Fig.  5). 


Fie.  5.  Kin.il  bind  structure  of  graphite  (tt  bands  are  shown  dashed). 
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6.  Discussion 

Most  of  the  previous  band  structure  calculations  obtain  very  similar 
values  for  certain  features  such  as  the  4.5  eV  for  M  point  n  transitions,  oc¬ 
cupied  a  bandwidth  of  nearly  16  eV,  and  it  -a  separation  of  around  12  eV. 
The  present  work  is  in  general  agreement  with  previous  theoretical  and  ex¬ 
perimental  work  as  displayed  in  Table  1.  In  particular,  the  gap  of  11.76  eV 
between  the  occupied  and  unoccupied  a  bands  at  T  and  the  separation  of 
11.3  eV  between  the  lower  a  and  it  bands  are  in  excellent  agreement  with 
the  experimental  results.  The  corresponding  values  of  the  muffin  tin  bands 
are  included  to  illustrate  the  dramatic  difference  between  these  and  final 
bands. 


Ackno  wledgmen  ts 

The  authors  thank  the  University  of  Pennsylvania  Graphite  Group  in 
general,  and  David  P.  DiVincenzo,  in  particular,  for  many  stimulating  dis¬ 
cussions  and  helpful  suggestions.  Special  thanks  are  due  to  the  staff  of  the 
Moore  School  Computing  Facility. 

This  work  was  supported,  in  part,  by  ARO  contract  DAAG-29^77-C400 
and  by  NSF  on  grant  number  DMR76-S0994. 


TABLE  1 

Comparison  of  graphite  band  structures 


KKR 

Final 

Willis 

Batra 

Zunger 

Exp. 

Ref. 

Hi] 

[10] 

[9] 

i.  rv-  rs* 

16  eV 

14. S  eV 

15.8  eV 

20.7  eV 

16.5 

gBW 

2.  To-  -  n* 

11.9 

11.76 

11.8 

12.5 

12.5 

11.5  i  1 

11 

occ.-unocc. 

3.  P3”  -  IV 

28 

18.63 

18 

15 

17.5 

77  BW 

4.  M3”  •  M4~ 

-4.4 

~4.3 

4.3 

4.5 

4.8 

4.6 

7 

M4  -M3” 

5.  T3*-  T2 

4 

1.65 

1.8 

1.3 

__ 

0.9  t  0.3 

22 

Trsep. 

6.  rr  -  T3* 

2.44 

11.3 

12 

16 

16 

11  »  0.3 

22 

7.  K2  -  K4 

3.2 

1.04 

1.3 

1.0 

1.6 

11 

Dw  “  bandwidth, 
scp.  ■  separation. 
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1.  INTRODUCTION 


The  electronic  transport  properties 
of  graphite  intercalation  compounds  have 
been  extensively  studied,1  with  emphasis 
on  the  metallic  conductivity  exhibited 
by  all  compounds  parallel  to  the  layers 
(the  a-direction) .  Alkali  metal  com¬ 
pounds  are  the  simplest  structurally, 
especially  the  stage  1  MCg  materials 
(M=K ,  Rb ,  Cs)  for  which  band  structures 
have  been  calculated. 7  Despite  their 
sensitivity  to  air,  these  compounds  are 
amenable  to  standard  galvanomagnetic 
techniques  so  that  one  can  obtain  reli¬ 
able  Hall  effect  and  magnetoresistance 
data  as  well  as  resistivity. 

Extensive  studies  on  well-character¬ 
ized  samples  based  on  HOPG  (highly 
oriented  pyrolytic  graphite)  were  first 
performed  by  Murray  and  Ubbelohde 3  (MU) . 
Their  interpretation  of  RH  data  for  MCg 
compounds  presumed  a  one-carrier  Fermi 
surface  which  was  subsequently  shown  to 
be  incorrect,4  and  their  analysis  of 
pa(T)  from  60-300K  in  terms  of  the  usual 
Bloch-Gruneisen  (BG)  formula  was  not 
convincing.  Suematsu  et.al.5(SHT)  mea¬ 
sured  Pa(T)  from  5-300K  and  found  a  good 
fit  to  0a (t)  =  A+BT+CT7,  and  proposed 
several  explanations  for  the  T?  term. 

Their  R„  data  conflicted  in  sign,  magni¬ 
tude  ana  temperature  dependence  with 
earlier  work  at  Penn,4'  and  their  re¬ 
sistance  ratios  R ( 300K) /R (v5K)  were 
considerably  lower.  In  a  series  of 
papers  the  Penn  group  reported  values  of 
pa»  and  magnetoresistance  at  several 
fixed  temperatures,4'6  along  with  exten¬ 
sive  measurements  of  transport  anomalies 
in  the  stage  2  MC 2g  compounds.6'7  Of 
these  3  sets  of  data,  the  MU  and  Penn 
results  agree  in  the  limited  range  where 
they  overlap,  but  the  SHT  and  Penn  results 
are  distinctly  different.  Thus  it  was 
deemed  advisable  to  repeat  some  of  the 
SHT  experiments. 

It  is  significant  that  both  MU  and 
the  Penn  group  used  i_n  situ  methods 
which  avoid  handling  the  samples  after 
compound  formation,  while  SHT  transfer 
their  samples  from  the  preparation  glass¬ 
ware  into  the  measurement  apparatus. 


The  only  quality  index  available  for 
comparing  samples  among  the  3  qroups 
is  R(300)/R(77)  for  KC„;  the  values  are 
11.4  (MU),3  5  (SHT)  b  and  11  (Penn).1* 

2 .  EXPERIMENTAL 

All  the  samples  studied  were  based 
on  HOP'S .  Rectangular  pieces  2.0  X  0.5  cm 
were  cut  using  a  diamond  impregnated 
wire  saw.  These  pieces  were  razor  cleaved 
to  a  thickness  of  approximately  0.5mm 
and  then  cut  into  five-arm  bridges  by  an 
air  abrasive  technique  to  minimize  inter¬ 
nal  damage.  The  five-arm  bridges  were 
then  tape  cleaved  to  a  thickness  of 
approximately  0.2mm.  The  distance  be¬ 
tween  voltage  arms  was  1.0cm  with  the 
Hall  arm  on  the  opposite  side  centered 
between  the  two  voltage  arms.  Electrical 
connections  were  made  with  platinum  wires 
and  sintered  gold  contacts.  The  samples 
were  then  mounted  in  a  specially  designed 
pyrex  rectangular  envelope  as  described 
previously. b  Vacuum-tight  seals  between 
the  Pt  wires  and  the  pyrex  were  achieved 
with  a  soft  intermediate  glass.  The  open 
end  of  the  rectangular  pyrex  was  then 
joined  to  a  long  pyrex  tube  into  which 
the  alkali  metal  could  be  distilled. 

After  distillation,  the  tube  containing 
the  alkali  metal  and  the  sample  was 
sealed  and  the  sample  intercalated  using 
the  usual  two-tempe rature  method.  Upon 
obtaining  the  desired  stage,  determined 
by  color  and  (00£)  x-ray  diffraction, 
the  sample  was  sealed  under  vacuum. 

Resistivity  measurements  were  made 
usinn  an  ac  current  of  10mA  at  33.3  Hz 
and  lock-in  detection  in  order  to  avoid  ther 
mal  emf  problems.  The  sample  temperature 
was  measured  with  a  calibrated  Si  diode 
and  controlled  with  an  electronic  temper¬ 
ature  controller  and  heater.  The  temper¬ 
ature  could  be  swept  up  or  down  at  a 
controlled  rate  automatically.  The 
voltage  droD  as  a  function  of  temperature 
was  recorded  by  a  semi-automated  system. 

Hall  data  were  obtained  by  sweeping  H  at 
fixed  t  . 


3. 


RESULTS 


The  “a  values  at  5K  and  300K  and 
their  ratios  are  listed  in  Table  I  for 
KCg ,  RbCg ,  and  CsGg.  For  KC  and  RbCg 
the  present  ratios  are  somewhat  better 
than  earlier  Penn  samples.  For  KCg  SHT 
report  one  sample  with  P^(l,2)  =  0.1 
li '1cm  and  three  samples  with  Pa(1.2)=0.1 
yQcm.  Thus  their  resistance  ratios 
range  from  3  to  60  times  less  than  ours, 
assuming  comparable  values  of  oa(300). 

The  temperature  dependence  of  Pa 
is  shown  in  Fig.  1  for  KC~.  Similar 
behavior  was  found  for  RbCg  and  CsCg, 
these  latter  two  being  essentially 
identical  and  slightly  more  resistive 
than  KC„ ,  as  was  found  by  MU.  It  is 
apparent  that  the  data  is  not  well- 
represented  by  the  usual  BG  formula, 
which  predicts  pVT5  at  low  T,  p'vT  at 
high  T.  Rather,  the  solid  line  in 
Fig.  1  is  a  linear  least  squares  fit  to 
ua (T)  =  A+BT+CT2  over  the  entire  temper¬ 
ature  range.  The  coefficients  are  list¬ 
ed  in  Table  II. 

In  contrast  to  the  stage  1  compounds, 
°a(T)  for  stage  2  compounds  is  quite 
linear 5~ 7 (except  for  the  anomalies)  over 
the  range  80-300K.  Here  again,  the  Penn 
samples  exhibit  lower  low- temperature 
resistivity  than  the  SHT  samples!  for 
KC^,  P  a  (1.2)  =2.4uf>cm  (SHT);  pa(4.2)  = 

0.25uficm  (Ref.  6).  This  may  be  signifi¬ 
cant  in  the  comparison  of  Hall  data  to 
follow . 

Figure  2  shows  the  field  dependence 
of  Rj,  (TT  |  |  c-)  at  low  temperature  for  KCg, 
CsCg  and  RbC^^.  It  is  clear  that  Rf,  is 
positive  and  varies  approximately  linear¬ 
ly  with  field  for  the  stage  1  compounds, 
while  R;}<0  and  is  field-independent  for 
RbC^^.  The  latter  is  consistent  with 
a  single-carrier  n-type  metal,  while  the 
former  suggests  more  complex  behavior. 

The  RbC.^  R^  value  is  independent  of  T 
up  to  300K,  while  RH  for  the  MC’s 

is  strongly  temperature-dependent,  as 
shown  in  Table  III.  The  data  newly  re¬ 
ported  here  are  consistent  with  earlier 
Penn  group  results4’6  in  showing  that 
for  MCg  R^  changes  sign  from  negative  at 


300K  to  positive  in  the  helium  ranqc , 
in  contrast  to  SHT.  The  samnle-to- 
sample  variations  among  the  three  Penn 
data  sets  are  smaller  than  the  varia¬ 
tions  observed  by  SI!T. 

Concerning  the  stage  2  compounds, 
Table  III,  both  groups  agree  that  Rjj 
is  negative  at  all  temperatures,  and 
that  the  temperature  dependence  of  the 
magnitude  of  Rjj  is  about  a  factor  of  2 
from  30 OK  to  ^4K  (KC24,CsC  ).  The 
Penn  values  are  however  consistently 
smaller  by  a  factor  of  3  or  4  than  the 
SHT  results.  The  magnitude  discrepancy 
between  SHT  and  Penn  results  gets  worse 
as  the  samples  get  more  dilute.  For 
stage  3  Onn  et.al.6  find  RH  = 

-21  X  10“4  cm3/ccul.  while  SHT  report 
-280  X  10-4;  the  respective  values  at 
stage  4  are  -57  X  10“4  and  either 
80  X  10“ 4  or  1000  X  10“4  (sign  not  deter¬ 
mined)  . 

4.  DISCUSSION 

The  ha  (T)  behavior  of  the  MCg 
compounds  is  clearly  unusual  for  such 
highly  conducting  materials.  It  is  ob¬ 
vious  from  the  Pa (T)  plots  that  no  T5 
or  linear-T  regions  occur  in  the  three 
stage-one  heavy-alkali-metal  GIC's  be¬ 
tween  5K  and  300K.  It  is  unlikely  that 
a  T s  region  can  be  found  below  5K  since 
the  residual  resistivity  seems  to  be 
dominant  even  at  5K.  It  would  not  be 
surprising  however  to  find  a  linear-T 
dependence  above  300K. 

The  reasons  for  the  inconsistency 
between  na(T)  and  the  BG  theory  are 
found  in  the  assumptions  which  lead  to 
the  BG  relation:  the  neglecting  of 

i)  separate  contributions  from  longitu¬ 
dinal  and  transverse  phonon  modes; 

ii)  electron-electron  scattering  and 

iii)  anisotropies  in  the  Fermi  surface. 
The  first  of  these  is  equivalent  to  ne¬ 
glecting  U-processes  which  are  dominated 
by  transverse  phonons.  General  carrier- 
carrier  scattering  cannot  be  overlooked, 
even  though  this  contribution  is  cnly 
significant  at  very  low  temperatures  for 
ideal  metals.9  Since  the  Fermi  surfaces 
are  partly  two-dimensional7  for  KCD  the 

O 


third  assumption  is  definitely  in  error. 

Table  II  indicates  the  temperature 
coefficients  resulting  from  fitting  to 
p  (T)  =  A  +  8T  +  CT;'.  The  first  term,  A 
represents  the  residual  resistivity  and 
8  is  related  to  electron-phonon  scatter¬ 
ing.  For  the  f  term,  two  possibilities 
must  be  considered.  First,  since  KCg 
has  partially  cylindrical  Fermi  surfaces, 
one  is  hopeful  that  Kukkonen's9  exten¬ 
sion  of  the  BG  theory  for  anisotropic 
metals  could  be  applied.  An  estimate  of 
when  the  T?  dependence  should  begin  can 
be  done.  For  N-processes  where  longitu¬ 
dinal  phonons  are  dominant,  the  onset 
temperature  for  T*’  behavior  is  T(T2)=230K 
Transverse  acoustic  phonons  cannot  be 
neglected  without  consideration,  however. 
The  transverse  sound  velocity  of  graphite 
( 2 70m/s )  yields  T(T?)  =  15K.  It  is  there¬ 
fore  important  to  find  out  the  relative 
weights  of  these  processes  in  the  com¬ 
pounds.  Kukkonen  did  not  have  this  pro¬ 
blem  with  bismuth  since  in  that  case 
TTA~TLA-  In  this  model  the  effective 
onset  temperature  must  be  far  less  than 
2  30K  for  KC„ .  The  cutoff  of  the  T? 
region  is  the  Debye  temperature,  234. 8K 
for  KCg.10  0D  is  however  an  average 
value ;°the  appropriate  characteristic 
temperature  may  be  above  300K.  The 
experimental  plots  show  no  purely 
linear-T  region  below  300K. 

In  light  of  these  facts,  and  since 
there  are  no  T 5  -+•  T?  -*■  T  transitions, 
a  better  explanation  of  the  T2  contribu¬ 
tion  might  be  the  second  possibility, 
general  carrier-carrier  scattering. 
Kukkonen  et  al.  [11]  also  presented  a 
theory  of  the  effect  of  collisions  be¬ 
tween  two  types  of  charge  carrier  on  the 
electrical  resistivity  of  a  simple  de¬ 
generate  semimetal.  Although  the  theory 
was  applied  to  TiS2  (a  layered  compound 
with  a  hexagonal  Brillouin  zone  and  aniso¬ 
tropy  of  10)  it  applies  in  general  to  any 
degenerate  electron-hole  system  with 
spherical  carrier  pockets  in  which  col¬ 
lisions  between  charge  carriers  conserve 
momentum.  The  result  of  this  theory  is 
a  T'  denendence  for  all  temperatures. 
Whichever  theory  is  more  applicable, 
modifications  will  have  to  be  made  to 


include  the  effects  of  the  two  differ¬ 
ent  electron  pockets  of  MCg  compounds. 

Several  observations  point  to 
carrier-carrier  scattering  rather  than 
anisotropy  as  the  origin  of  the  T: 
behavior  in  MCg  compounds.  First,  the 
T2  contribution  is  much  stronger  in  stage 
1  than  in  stage  2  compounds  (the  effec¬ 
tive  exponents  of  log-log  plots  are 
respectively  1.8  and  ^1  over  the  range 
60-300F) ,  whereas  stage  2  is  more  highly 
anisotropic.1  Second,  the  magnitude  of 
the  T2  coefficient  decreases  in  the 
sequence  K:Rb:Cs,  as  does  the  density  of 
electrons  in  the  quasi-spherical  pocket,12 
suggesting  that  the  T2  term  might  origi¬ 
nate  in  collision  between  electrons  from 
the  spherical  pocket  and  electrons  (or 
holes)  originating  from  the  cylindrical 
Fermi  surface  pockets. 

The  8  and  C  coefficients  obtained 
for  KC„  are  in  general  agreement  with 
the  SHT  results,  although  they  had  some 
samples  with  much  larger  linear-T 
terms.  The  major  discrepancy  lies  with 
the  defect  term  A,  as  discussed  previous¬ 
ly- 

A  major  effort  was  put  forth  in 
this  work  to  obtain  values  for  stage 
1  and  2  compounds,  in  light  of  the 
discrepancies  between  the  recent  SHT 
data5  and  earlier  Penn  results.4  Once 
again  we  found  R^  >0  at  low  T  for 
stage  1  compounds,  in  contrast  to  SHT. 

The  discrepancy  may  be  due  to  a  pre¬ 
dominance  of  defect-related  phenomena 
in  the  SHT  samples,  which  in  turn  may 
be  related  to  the  difference  between 
in  situ  and  erx  situ  techniques. 

The  sign  change  from  ^0  at  high  T 
to  0  at  low  T  suggests  that  holes 
dominate  at  low  T.  The  calculated  Fermi 
surface  does  contain  hole  orbits  through¬ 
out  an  appreciable  portion  of  the  zone, 
due  to  the  connectivity  of  quasi-spheri¬ 
cal  and  cylindrical  electron  pockets.2 
The  hole  orbits,  centered  at  M  and  Pi 
points,  show  very  small  curvature  over 
part  of  the  trajectory,  and  thus  may 
have  very  high  mobility  at  low  T. 

Another  possible  reason  for  the  cross¬ 
over  from  electron  to  hole  dominance 


could  be  a  negative  thermal  expansion 
coefficient  (as  for  graphite  itself), 
in  which  case  the  basal  plane  contraction 
of  the  zone  upon  cooling  would  increase 
the  fraction  of  c*  over  which  the  two 
types  of  electron  pockets  are  connected. 
This  in  turn  would  enhance  the  fraction 
of  hole-type  orbits.  Note  that  the  co¬ 
existence  of  electron  and  hole  orbits 
evidenced  by  would  be  consistent  with 
electron-hole  scattering  as  the  origin 
of  the  T'  term  in 

A  second  point  of  contention  between 
the  Penn  and  SHT  results  is  whether  or 
not  for  stage  2  compounds  is  consistent 
with  a  single-carrier  model.1'4  The  metal 
atom  densities  are  ^3,5  X  1021/cm3  in 
MC24.  (M=K,  Rb ,  Cs)  ,  so  the  minimum 
absolute  value  of  Rh  (equivalent  to  one 
free  electron  per  M  atom)  is  ^20  X  10~4 
cm3/coul.  Referring  to  Table  III,  this 
limit  appears  to  be  satisfied  by  KC24 
and  CSC24  at  4K,  but  not  by  RbC24, 
at  any  T  nor  by  the  others  at  T>77K. 

A  slight  departure  from  ideal  one-carrier 
behavior  might  be  simply  due  to  the 
"trigonal  warping"  of  the  cylindrical 
electron  surfaces.2  Of  the  three  com¬ 
pounds,  only  RbC24  exhibits  an  RH  which 
decreases  in  magnitude  with  decreasing 
T,  resulting  in  a  major  departure  from 
the  one-carrier  limit  at  4K. 

The  differences  between  electron/ 
atom  ratios  (or  fractional  ionization) 
inferred  from  the  Penn  and  SHT  data  are 
substantial.  For  KC24  f=l.l  (Penn)  or 
0.3  (SHT);  for  KC36  f=l.l  (Penn)  or 
0.1  (SHT);  for  KC4g  f =0 . 6  (Penn)  or  0.03- 
0.3  (SHT).  We  also  find  f^l  for  CSC246 
and  CsC3g.4  Many  other  experiments  are 
consistent  with  f=l  for  stage  2  and 
higher  alkali  metal  compounds,1  We 
suggest  again  that  the  difference  may 
arise  from  differences  in  sample  handling 
technique . 

Since  the  galvanomagnetic  properties 
are  so  complex,  their  detailed  interpre¬ 
tation  will  have  to  incorporate  the  sub¬ 
tle  features  of  the  Fermi  surfaces. 
Establishing  the  latter  requires  testing 
the  band  structure-derived  Fermi  surface 


with  quantum  oscillatory  experiments 
such  as  the  de  Haas-van  Alphen  effect. 

SHT  have  performed  such  experiments  on 
the  same  specimens  used  in  their 
galvanomagnetic  studies.13  The  values 
of  f  inferred  from  the  dHvA  results  are 
M) .  2  for  KC^g  and  KC^g.  Similar  ill  situ 
experiments  on  our  samples  are  in  progress 
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TABLE  II 


Basal  resistivity  ‘ a  at  300K  and  5  K 
for  stage  1  alkali  metal  compounds. 

The  ratio  is  a  rough  index  of  sample 
quality.  Ratios  from  previous 
publications  are  included. 

Parameters  obtained  from  fitting 
Pa  =  A  +  8T  +  CT?.  Data  from  Ref. 

14  is  included  for  KCg. 

Hall  coefficients  measured  at  ^20  kO 
at  several  temperatures  for  a  variety 
of  alkali  metal  compounds.  Earlier 
data  are  included  for  comparison. 
Values  in  parentheses  are  upper  limits 
only.  The  sign  of  R^  for  KC^g  was  not 
determined  in  Ref.  5; 


Figure  Captions 
Figure  1 


Figure  2 


Basal  resistivity  Pa  versus  T  for 
KCg  (i_n  situ  4-wire  measurement)  . 
The  solid  line  is  a  linear  least 
squares  fit  to  Pa  -  A  +  fiT  +  CT2 
(fit  parameters  listed  in  Table  IT) 
The  semilog  scale  exaggerates  the 
discrepancy  between  data  points 
and  the  linear  fit  below  50K. 

Hall  coef_f icient_versus  magnetic 
field  (  3"  |  |  a,  H  |  |  c  )  for  KC^, 
CsC„,  and  at  low  temperature. 

With  a  constant  voltage  error,  the 
RH  error  decreases  with  increasing 

H. 
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HALL  COEFFICIENT  (10  cm 


